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1 ■ ABSTRACT 

o 

(NT Context. The unified scenario for Active Galactic Nuclei (AGN) postulates that our orientation with respect to a parsec- 

scale azimuthally-symmetric gas and dust system causes their different phenomenology in the optical/UV and X-ray 

■ bands. Only recently, high-resolution radio (VLBI) and IR interferometric observations have provided direct constraints 
' on size and structure of this obscuring system (known historically as the "torus"). On the other hand, variability of 

optically-thick X-ray absorption and reprocessing in heavily obscured AGN often probe smaller scales, down to the 
Broad Line Region and beyond. 

Aims. We aim at constraining the geometry of the reprocessing matter in the nearby prototypical Seyfert 2 Galaxy 
Markarian 3 by studying the time evolution of spectral components associated to the primary AGN emission and to its 
Compton-scattering. 

Methods. We have analyzed archival spectroscopic observations of Markarian 3 taken over the last ~12 years with the 
XMM-Newton, Suzaku and Swift observatories, as well as data taken during a monitoring campaign activated by us in 
^ ■ 2012. 

Results. The timescale of the Compton-reflection component variability (originally discovered by ASCA in the mid- 
'90s) is jS64 days. This upper limit improves by more than a factor of 15 previous estimates of the Compton-reflection 
variability timescale for this source. When the light curve of the Compton-reflection continuum in the 4-5 keV band 
is correlated with the 14-195 keV Swift/BAT curve a delay ^1200 days is found. These results are dependent on the 
model used to fit the spectra, although the detection of the Compton-reflection component variability is independent 
of the range of models employed to fit the data. Reanalysis of an archival Chandra image of Markarian 3 suggests that 
the Compton-reflection and the Fe K a emitting regions are extended to the North up to ~300 pc. The combination of 
these findings suggests that the optically-thick reprocessor in Markarian 3 is clumpy. 

Conclusions. There is mounting experimental evidence for the structure of the optically-thick gas and dust in the nuclear 
\Q , environment of nearby heavily obscured AGN to be extended and complex. We discuss possible modifications to the 

standard unification scenarios encompassing this complexity. Markarian 3, exhibiting X-ray absorption and reprocessing 

■ on widely different spatial scales, is an ideal laboratory to test these models. 
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OV 

o: 

^ | 1. Optically thick reprocessing in AGN This evidence suggests that one is dealing with a structure 

• • ■ at a parsec-scale distance from the active nucleus. 

■ Unified scenarios for radio-quiet Active Galactic Nuclei 
(AGN) ()Antonucci k Miller 19851 lAntonucci 1993D postu- The unlfied scenarios have been shown to hold for 
late the existence of an azimuthally-symmetric distribu- asma 2L _as_97% of AGN in X-ray selected samples 

& ■ tion of optically-thick gas and dust surrounding the nu- (|Mateos et al. 2 010J). Violations to its basic prediction (i.e.: 

1 cleus, as well as the gas clouds responsible for optical broad a different classification from optical and X-ray spec- 

(£1000 km/s) lines. This structure is commonly referred as tr0SC0 Py) c f n be explained by the poor quality of the spec- 

the "torus" . We will adhere to this convention in this paper, trum 4 011 which the classification was based, variability of 

as customary in the recent AGN literature, although con- the AGN emission complexity m the absorber structure 

straining the geometry of this absorbing system is one of the (|Bianchi et al. 2 .012J , or threshold effects on the AGN ly> 

main astrophysical motivation of the study presented in this mmosity. For instance, a minimum Eddmgton ratio could 

paper. Warm dust dominates the mid-IR emission in nearby be necessary fo r the Broad Line Regions (BLR) to form 

AGN (see, e.g. jAlonso-Herrero et al., 20111 The innermost dNicastro 2000[ | Marmucci et al. 2012a) 
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regions of the dusty phase has to be therefore located be- It is crucial to validate the above structure model by 
yond the sublimation radius (~ 0AL 1 / 2 pc, where L 45 is th e constraining the geometry of the gas and dust in the AGN 
ionizing luminosity in units of 10 45 erg s'^ IBarvainis 1987|) . environment. Reverberation ma pping of optical- UV lines 

against their ionizing continuum (Pete rson et al. 2004P and 
monitoring of changes in the column density covering 
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therein) constrain the location of the BLR clouds within a 
few light-days in typical Seyfert-like objects (i.e., with total 
bolometric luminosity < 10 44 erg s _1 ). However, it is cur- 
rently arguable whether X-ray obscuration can be indeed 
ascribed to BLR clouds in most obscured AGN, in partic- 
ular in Compton-thick ones {N H > <r t -1 =1.6 x 10 24 cm- 1 ). 

Recently, mid-IR interferometry has allowed for the 
first time to effectively "image" the dusty component 
of the parsec-scale torus. In the three cases for which 
such measurements have been possible so far (NGC 1068, 
the Circinus Galaxy, and NGC 4151), the data are 
consistent with a two-temperature distribution of dust, 
with most the of the obscuration occurring on scales 
smaller than 1 pc jJa f fe et al. 20041 |Tristram et al., 2007[ 
IBurtscher et al. 2009p . 

Constraining the geometry of the X-ray reprocessing gas 
is currently possible only by studying the response of spec- 
tral components dominated by optically-thick reprocessing 
to variability of the primary continuum. These measure- 
ments are difficult, because they require simultaneous mea- 
surements of fluxes in at least two energy bands dominated 
by the AGN primary emission and by its reprocessing, re- 
spectively, on times scales long enough to be able to de- 
tect the response of the latter to changes of the former. 
These timescales are likely to be of the order of months, at 
least. Dedicated monitoring campaigns of this kind are ex- 
pensive, although the launch of the Swift satellite in 2004 
(jGehrels et al. 2004[) . with its BAT X-ray transient mon- 
itoring instrument scanning continuously the high-energy 
sky in the 15-150 keV band, has opened a new window in 
this field. 

In this paper we discuss flux variability in a Compton- 
reflection dominated energy band on timescales of months- 
years in the nearby (z=0.014) prototypical Seyfert 2 galaxy 
Markarian 3. Known since the dawn of X-ray astronomy, 
Markarian 3 hosts a heavily obscured AGN. The most 
prominent feature in its X-ray spectrum is a strong iron 
K a fluorescent line ([Awaki et a l. 1991|) . The continuum in 
the 3-7 keV range is dominated by unabsorbed Compton- 
reffection ( Gcorgantop oulos et al. 19 99, Cap pi et al. 1999[ ). 
Iwasawa et al. (1994) discovered a variability of the hard 
(E >4 keV) luminosity by a factor of 3 in 3.6 years by com- 
paring their ASCA observation in April 1993 with previous 
Ginga and BBXRT observations. The soft X-ray emission 
had remained constant over the previous 13 years within 
statistics. 

The paper is organized as follows: in Sect. 2 we present 
the observations discussed in this paper, as well as their 
data reduction. The spectral analysis is described in Sect. 3, 
where we also show light curves over the first decade of this 
century in energy bands dominated by different spectral 
components. We discuss the implication of our results in 
Sect. 4, focusing on three main aspects: the geometry of 
optically-thick reprocessing matter; the location of the iron 
K„ line emitting gas; and the variability of the soft X-ray 
emission. We summarise our main results in Sect. 5. 

In this paper: energies are shown in the source's rest 
reference frame; statistical errors are calculated at the ler 
confidence level for one interesting parameter; and the fol- 
lowing cosmological parameters were used to calculate lu- 
minosities: H o =70 km s" 1 Mpc" 1 , A o =0.73, and O A /=0.27 
(jBennett et al. 2003p . At the distance of Markarian 3, 1" 
corresponds to 270 pc. 



2. Data reduction 

In this paper we analyse Markarian 3 spectra obtained with 
imaging CCD cameras on-board Swift and XMM-Newton. 
The longest XMM-Newton obser vatio n (October 2000) is 
discussed in Bia nchi et al. 20051 and IPounds et al. 2005. 
The log of the observations discussed in this paper is shown 
in Tab. [Q 

2.1. XMM-Newton /EPIC-pn 

In this paper we discuss only spectra extracted with the 
EPIC-pn camera (Str iider et al. 20"0lj) . Markarian 3 was ob- 
served seven times in the early phase of the XMM-Newton 
mission using the Full Frame Mode. The March 2012 obser- 
vation was instead performed in Timing Mode. Raw scien- 
tific telemetry (Observation Data Files) was reduced using 
SAS Version 12 (Gab riel et al. 2003)) . and the calibration 
files available on June 1 2012. Spectra were extracted from 
calibrated event lists generated with the task epchain. The 
Timing Mode event list energy scale was corrected for any 
rate-dependent effect, which is, however, likely to be negli- 
gible due to the source flux level. 

For imaging modes, source spectra were extracted 
using circular regions of 40" radius centered around 
Markarian 3 optical coordinates (a20oo=93. 9015194°; 
#2000=71.0375254°); background spectra were extracted 
from circular regions free of serendipitous contaminating 
sources on the same CCD quadrant as the boresight posi- 
tion, and corresponding to RAWY positions close to where 
Markarian 3 is located on the EPIC-pn boresight CCD, to 
ensure that a similar Charge Transfer Efficiency correction 
applies. 

For the Timing Mode observation, source spectra were 
extracted from a box in detector coordinates centered on 
RAWX=36 and with a radius of 5 pixels (each pixel in the 
EPIC-pn camera corresponds to 4.4"). No constraints were 
applied to the pseudo-spatial coordinate RAWY, which rep- 
resents the elapsed time during the observation in this 
mode. Background spectra were extracted from a nearby 
box centered on RAWX=50 to ensure that the background 
spectrum represent as faithfully as possible the true back- 
ground underneath the source position. The background 
represents about 75% of the source level in the 4-10 keV en- 
ergy band. The emission profile along the RAWX coordinate 
shows a gradual decrease (once the signal corresponding to 
Markarian 3 is removed) by about 5% between RAWX=27 
and RAWX=55. Through a linear interpolation to RAWX=36, 
we estimate that our background spectrum could underes- 
timate the true background level underneath the source by 
about 3%. We have added this factor in the analysis of the 
Timing Mode spectrum. 

Response files appropriate for each observation were 
generated with the SAS tasks rmfgen and arfgen. The 
EPIC-pn response is known to be stable within 3% 
(jSartore et al. 20121) . 

2.2. Swift/XRT 

Swift/XRT ([Burrows et al. 2005[) observed Markarian 3 on 
9 occasions. However, as 2 observations have relatively short 
exposure time (2 ks and 55 s respectively), we consider here 
7 Swift XRT observations with exposure times longer than 
5 ks. We extracted Swift XRT products from the cleaned 
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Table 1. The log of the observations analyzed in this paper. The nomenclature "[nm]" indicates the Swift/XKT obser- 
vations which have been merged together. For instance: "0003546000 [23]" means that a single spectrum was extracted 
from Obs.#00035460002 and Obs.#00035460003. 



Obs.# 


Start Time 


(ks) 


Instrument 


Mode 


0111220201 


2000-10-19T15:39:46 


37.0 


EPN 


PrimcFullWindow 


0009220601 


2001-03-20T20:49:05 


4.9 


EPN 


PrimeFullWindow 


0009220701 


2001-03-28T21:33:43 


3.0 


EPN 


PrimcFullWindow 


0009220901 


200 1-09- 12T01: 12:44 


1.7 


EPN 


PrimeFullWindow 


0009221601 


2002-09-16T05:32:03 


1.6 


EPN 


PrimeFullWindow 


0009220401 


2002-03-10T13:57:37 


2.6 


EPN 


PrimeFullWindow 


0009220501 


2002-03-25T17:45:25 


3.4 


EPN 


PrimeFullWindow 


00035460001 


2006-03-21T07:16:01 


5.1 


XRT 


PHOTON 


0003546000 [23] 


2006-04-02T12:00:00 


18.1 


XRT 


PHOTON 


00037226001 


2008-01-22T00:40:08 


9.4 


XRT 


PHOTON 


0003546000 [45] 


2012-01-11T12:00:00 


14.7 


XRT 


PHOTON 


0656580301 


2012-03-15T12:32:16 


30.9 


EPN 


FastTiming 


00035460008 


2012-05-06T03:18:12 


9.2 


XRT 


PHOTON 



photon counting event files following the procedures de- 
scribed in The SWIFT XRT Data reduction Guide V.lM 
Spectra of Markarian 3 were extracted from circular re- 
gions of radius 60", while background spectra were obtained 
from larger (100" in radius) source-free regions close to the 
source. We used the standard response files for the photon 
counting mode with grades to 12 which are part of the cur- 
rent Swift XRT calibration database, taking into account 
the observation epoch. Ancillary files have instead be cre- 
ated for each observation using the xrtmkarf task for the 
specific source position. After checking that the individual 
spectra were consistent, we have merged the spectra from 
ObsID 00035460002 and 00035460003 (3 days apart) and 
from ObsID 00035460004 and 00035460005 (1 day apart). 
Hence 5 X-ray XRT spectra of of Markarian 3 are consid- 
ered here. 



2.3. Swift/ BAT 



The Swift-BAT ( |Barthelmy et al. 2005"]) data were pro- 
cessed with the Bat .Imager code ( Segreto et al. 2010 ). The 
15-150 keV light curve, as well as those in three energy 
bands: 15-30, 30-70, and 70-150 keV, were obtained pro- 
ducing a set of all-sky maps each covering 15 days, starting 
from December 2005 up to March 2012, and extracting the 
count rate and its error from the pixel corresponding to the 
source position in each of the maps. 



3. Spectral analysis 

All the EPIC-pn and XRT spectra were independently 
analysed in the nominal 4-10 keV band, where Compton- 
reflection dominates fllwasawa et al. 1994|) . However, the 
upper bound is observation-dependent (see Fig. [T0|) as XRT 
dost not confidently detect the source up to 10 keV in 
most observations. Background-su btracted sp ectra were fit- 
ted using the Cash statistics (C; ICash 19761) f° r us to be 
able to use the natural instrumental spectral binning. 



http: / / swift .gsfc.nasa.gov /docs / swift /analysis / 



Our first goal was to measure the light curves of the 
fluxes of spectral components associated with reflection of 
the primary continuum by optically thick matter. We have 
therefore fit the spectra with the following phenomenolog- 
ical model, M pr ( "pure reflection" ) : 

M pr = C R (T,E) + ^ 1 G i (E) 

where Cr is a Compton-reflection continuum component 
(Magdziarz & Zdziarski 1995) and the Gi are unresolved 
(i.e., intrinsic width assumed to be equal to 0) Gaussian 
emission lines. In the Compton-reflection component we 
have assumed: no high-energy cut-off of the primary con- 
tinuum (as this will most likely fall outside the EPIC-pn 
and XRT sensitive bandpass: Dadma 2l)l)7i: solar abun- 
dances (according to Anders & Grevesse 1983); and a 45° 
inclination angle between the normal to the plane of the 
Compton-reflecting slab and the line-of-sight. The only free 
parameter in this component, apart from its normalization, 
is therefore the photon index of the power-law primary con- 
tinuum, r. The centroid energies of the Gaussian profiles 
were fixed to 6.4, 6.7, and 6.96 keV, respectively, corre- 
sponding to fluorescent emission from neutral or mildly 
ionized iron, the resonant component of the He-like triplet 
and the H-like iron line, respectively. Sako et al. (2000) 
report for the Hea an energy of 6.685 keV (without er- 
ror) from their analysis of the Chandra /HETG spectrum of 
Markarian 3. This energy is intermediate between that of 
the resonant, the intercombination (6.668-6.682 keV) and 
the forbidden (6.637 keV) component. Neither the EPIC- 
pn nor the Swift/XRT spectra can constrain this param- 
eter. We have therefore opted for fixing the centroid en- 
ergy to the laboratory value of the resonant component 
following the results after Bianchi et al. (2005), which rule 
out an important contribution of the forbidden component. 
The results presented in this paper are insensitive to this 
choice. The intensity of the He-like and H-like lines can 
be constrained only in the longest XMM-Newton observa- 
tion: / 6 . 7 =(7.0±1.4)xl0- 6 and 7 6 . 96 =(5.5±1.4) x 10~ 6 pho- 
tons cm~ 2 s , respectively. In the other observations they 
have been fixed to the above best-fit values. This assump- 
tion is astrophysically justified, because a substantial frac- 
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Table 2. Best-fit values of the free parameters in M pr (T, £4-5, Iq.4, C/v) and of the phenomenological model used to 
fit the 0.5-2 keV energy band (T s , Fo.5_i.g). -F4-5 (-F0.5— 1.5) is the flux of the Compton-reflection component (total flux) 
in the 4-5 (0.5-1.5) keV energy range; Iqa is the normalization of the iron K a line; C is the best-value Cash statistic 
value. 





r 




-Fo.5-1.5 a 


F 4 -5 a 


T b 
16.4 


C/v 


0111220201 


1.00 ± 0.04 


2.79 ±0.03 


-12.112 ±0.004 


-12.38 ±0.01 


4.07 ±0.17 


1338.6/1289 


0009220601 


1.19 ± 0.014 


2.85 ± 0.08 


— 12.117 ± 0.011 


— 12.48 ± 0.04 


4.0 ± 0.5 


575.9/769 


0009220701 


1.2 ± 0.2 


2.78 ±0.11 


-12.109 ±0.014 


-12.55 ±0.06 


4.2 ± 0.6 


367.0/533 


0009220901 


1.1 ±0.3 


3.14 ±0.15 


-12.110 ±0.019 


-12.53 ±0.08 


3.2 ±0.8 


248.3/368 


0009221601 


1.4 ±0.3 


3.17 ±0.15 


-12.090 ±0.019 


-12.58 ±0.08 


3.9 ± 1.0 


239.7/326 


0009220401 


1.3 ± 0.2 


2.89 ±0.12 


-12.121 ±0.015 


-12.51 ±0.06 


3.4 ±0.7 


404.0/520 


0009220501 


1.5 ±0.2 


2.99 ±0.10 


-12.100 ±0.013 


-12.52 ±0.05 


3.3 ±0.5 


409.2/580 


00035460001 


1.1 ±0.8 


2.8 ±0.4 


-12.18 ±0.05 


-12.42 ±0.11 


3.4 ± 1.7 


70.5/86 


0003546000 [23] 


1.7 ±0.3 


2.60 ±0.16 


-12.14 ±0.02 


-12.37 ±0.05 


2.5 ±0.8 


213.6/259 


00037226001 


0.8 ±0.8 


2.8 ±0.3 


-12.25 ±0.04 


-12.66 ±0.11 


3.1 ± 1.2 


78.4/102 


0003546000 [45] 


1.4 ±0.3 


2.76 ±0.14 


-12.16 ±0.03 


-12.37 ±0.06 


5.1 ± 1.7 


202.2/242 


0656580301 


0.8 ±0.3 


3.02 ±0.10 


-12.163 ±0.010 


-12.54 ±0.08 


4.0 ±0.5 


1314.8/1296 


00035460008 


1.4 ±0.5 


2.4 ±0.2 


-12.25 ±0.03 


-12.53 ±0.09 


3.3 ± 1.3 


95.2/127 



a base-10 logarithm of the flux in units of erg cm 2 s 
6 in units of 10 5 photons cm -2 s _1 



tion of them is likely to produced in the X-ray Narrow Line 
Regions, extended on scales of the order of a few hundred 
parsecs (see Bianchi et al. 2006 for a quantitative discussion 
of this point). M pr does not leave unaccounted for residuals 
in any of the spectrum. The best-fit parameter values are 
reported in. Tab. [2] 

Fig. [1] (upper and medium panels) shows the light curves 
of the observed flux in the 4-5 keV band, and of the inten- 
sity of the K a fluorescent iron line. We have also added the 
results of the Suzaku pointed observation of Markarian 3 
discussed by Awaki et al. (2008; Obs.#1000040010)Q We 
decided not to analyze again the XIS spectra, because 
they were reduced and published by a highly experienced 
team, and no significant change in the calibration occurred 
above 3 keV since the time of the original data reduction 
in Awaki et al. (2008). The 4-10 keV energy band, cho- 
sen because it is dominated by emission-line free Compton- 
reflection, shows a variability dynamical range larger then 
70%. Similar dynamical ranges are seen when comparing 
the XMM-Newton or Swift observations among each other, 
ruling out cross-calibration effects (which should be anyhow 
of the order of at most 10%; |Tsujimoto et al. 20lT| ). The 
chi-squared for a fit against a constant, Xcl v i is 199.9/13. 
On the other hand, the iron line light curve exhibits a lower 
variability dynamical range. The Xc/ U value is 34.1/13. The 
two quantities are very weakly correlated (Fig. [5]). A lin- 
ear fit I /{I) = a + b x F/(F) yiekkfl a = 0.7 ± 0.2 and 



2 Markarian 3 was observed also by the CTiandra/HETG in 
combination with ACIS-S exactly seven months before the first 
XMM-Newton observation. We downloaded the spectrum from 
the TGCat archive, and analyzed it using the same procedure 
and model as the EPIC-pn and XRT spectra. While the in- 
tensity of the iron K a line is comparatively well constrained 
[J 6 .4=(5.0 ± 0.5) x 10" 5 erg cm 2 s'^ ISako et al. 2000| . the 4- 
5keV Compton flux is not (log 10 (F 4 _ 5 )=-12.8±0.7; cf. Tab.©. 
We will therefore not discuss these data in our paper. 

3 We normalize here I and F to their mean in order to reduce 
the errors on the best-fit parameters, which otherwise could be 



Mkn3 



r "'^^f^i^ 









1X10 13 2X10" 13 3x10" 13 4x10' 13 5X10" 13 6x10 13 

4-5 keV flux (cgs) 

Fig. 2. Intensity of the Fe K a against the observed flux in 
the 4-5 keV energy band. The dotted line indicates the locus 
of constant line EW, normalized to the value measured dur- 
ing the deepest XMM-Newton observation (October 2000). 
The dashed line indicates the linear best-fit, whereas the 
dot-dashed lines indicate the envelope corresponding to the 
la uncertainties on the best-fit parameters. 

b = 0.32±0.17. The two-side significance of the Spearman's 
rank correlation coefficient is ~31% only. This indicates 
that Equivalent Width (EW) is not constant. However, care 
must be exercised in interpreting this result, as well as the 
low x\/v of the Fe K Q intensity light curve, because the 
statistical uncertainties on its measurements are a factor 
of 2-3 larger than those on the 4-5 keV Compton flux. To 
make a more quantitative assessment of the correlation like- 
lihood (or lack thereof) between the observables in Fig. [2l 
we performed Monte-Carlo simulations of the Fe K Q in- 
tensity light curve, assuming that it follows the same rela- 



artificially increased by the leverage induced by the large differ- 
ence in absolute value between the fitted quantities 
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Fig. 1. Light curves of the: observed flux in the 4-5 keV energy band (upper panel); intensity of the Fe K a line (medium 
panel); observed flux in the 0.5-1.5 keV energy band (lower panel). The numbers in each panel indicate the chi-squared 
value for a ft against a constant versus the number of degrees of freedom. Filled circles: XMM-Newton/EPIC; empty 
circles: Swift/XRT; filled star. Suzaku. The scale on the y-axis in all panels corresponds to ±50% of the mean level 
(dotted lines). 



tive trend (in units of ratio against the mean) as the mea- 
surements of the 4-5 keV flux. In astrophysical terms, the 
simulated scenario assumes that the continuum Compton- 
reflection and the Fe line intensity simultaneously respond 
to changes of the same primary continuum. For each ob- 
servation, 10 4 realizations of a Gaussian distribution where 
drawn from a random sample, whose mean is the Fe line 
intensity aligned to the Compton-flux, and whose standard 
deviation is the measured statistical uncertainty on the Fe 
line intensity. For each simulated light curve, we calculated 
the slope of the best linear fit between the Compton flux 
and the simulated Fe line intensity. In Fig. [3] we compare 
the distribution of the simulated slopes with the measured 
slope. The probability of getting a value of b < b is <10~ 4 . 

Variation of the flux in the 4-5 keV energy band could 
be due to a change of the Compton-reflected flux in this 
band, as well as to the emerging of different spectral com- 
ponents. It is in principle possible, for instance, that a re- 
duction of the photoelectric column density covering the 
primary AGN emission brings a fraction of the primary 
emission into this energy band during some of the obser- 
vations presented in Fig. [T] The intrinsic spectral indices 
yielded by M pr are significantly flatter than typically ob- 
served in Seyfert galaxies (jBianchi et al. 2009aP and mea- 
sured in X-ray broad-band observations of Markarian 3 
( Cappi et al. 1999| lAwaki et al. 2008|l . While we warn 
against interpreting T in M pr as a reliable measurement 
of the photon index of the primary continuum on spectra 
with a small sensitive bandpass, the systematic flat spec- 
tral indices may indicate that a different spectral model is 
required to properly fit the data. 




Slope 

Fig. 3. Distribution of slopes of the best liner fit between 
the 4-5 keV Compton fluxes and Fe line intensities, when 
the latter are simulated assumption that they responded 
simultaneously to the same primary emission. The vertical 
line indicates the measured slope b. 

In order to test the impact that different origins for 
the spectral variability may have on our measurements of 
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the Compton-reflection variability, we fit the EPIC-pn and 
XRT spectra with two alternative models. These models 
have been applied to the EPIC-pn and XRT spectra dis- 
cussed in this paper under the same conditions as M pr . 
In particular, we fit the data in the same energy band: 4— 
10 keV. Because we want to address the flat spectral indices 
yielded by M pr , in the reminder of this Section we will fix 
the value of this parameter to that determined by Awaki 
et al. (2008; Ti nt = 1.8) from their analysis of the Suzaku 
observation of Markarian 3, unless otherwise specified. This 
choice is justified by the fact that the spectral coverage of 
the scientific payload on board Suzaku up to several tens 
of keV allows one to better constrain the intrinsic shape of 
the primary continuum in such an obscured source. 

- unobscured reflection plus obscured transmission. We as- 
sumed the full continuum model as discussed in 
IBianchi et al. 20051 

M urot = C R {T mt ,E) + Y? l=1 G l {E) + e-°^ NH -r x AE T ^ 

There is a strong correlation between the primary spec- 
tral index and the iron abundance of the Compton- 
reflection component (see, e.g., the discussion in Bianchi 
et al. 2005), primarily driven by the depth of the neu- 
tral iron photo-absorption edge. We therefore imposed Zp e 
to vary within the confidence interval as determined from 
the deepest XMM-Newton observation: Z Fe =O.43±O.O9Z . 
The photoelectric cross-section a p h(E) is according to 
the Tu bingen-Boulder mo del with the associated abun- 
dances (Wilms et al. 2000), and includes also the contri- 
bution from optically-thin Compton-scatterin^ The typ- 
ical column densities measured in the EPIC-pn spec- 
tra are iV ff ,j,~1.2-1.6xl0 24 cm" 2 . In the deepest EPIC- 
pn observations:iVH,p=(1.19 ± 0.06), and (0.49±g;Jg) x 
10 24 cm" 2 in October 2000 and March 2012, respectively. 
Lower limits from the XRT spectra are in the range 
_/V# )P ~0.6-1.4x 10 24 cm -2 . We urge readers not to over- 
interpret the measurements of the absorber column density 
in this scenario. The column density of the absorber cannot 
be constrained by our data if the value of the photon index 
is left free to vary. Fixing the shape of the spectral index 
severely limits the size of the confidence interval parameter 
space. 

- obscured reflection . We added a photoelectric absorption 
component covering the Compton-reflection emission: 

M or = e-^"^ N ^C R (T int ,E) + ^GiiE) 

We again constrained Zi? e =0.43±0.09Z Q . In astrophysi- 
cal terms, M or represents a scenario whereby the line-of- 
sight to the inner far side of the torus grazes the rim of 
its near side or crosses a dust lane in the host galaxy 
(|Lamastra et al. 2008[) . Although gas from the torus rim 
is likely to be evaporated by the AGN radiation pressure, 
and therefore significantly ionized, we used in M or a pho- 
toelectric model for cold matter, as the data do not require 
the ionization parameter as a further degree of freedom. 

4 We employed in this context the implementation cabs in 
Xspec for comparison with previous works, despite the issues 
described in Murphy & Yaqoob (2009) which should discourage 
its use for Compton-thick absorbers. They primarily affect the 
determination of the absorption-corrected luminosity of the ob- 
scured primary, a measurements in which we are not interested 
for the scope of this paper. 
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Energy (keV) 

Fig. 4. Spectrum of the XMM-Newton October 2000 
observations (upper panel) and residuals in units of 
data/model ratio (upper panel) if a model of ionized re- 
flection with rj„t=1.8 is applied. 



We refrain from deriving any astrophysical inference from 
the measured column densities, due to the simplicity of the 
assumptions and to the statistical quality of most of the 
spectra analyzed in this paper. Not surprisingly, only the 
two deepest XMM-Newton observations require a column 
density significantly different from zero: JV/, jr =4.5±0.6 and 
(17±|) x 10 22 cm" 2 in October 2000, and March 2012, re- 
spectively. In the other observations, the statistical quality 
of the spectra is the main driver behind upper limits in the 
range 1.7-6 xlO 22 cm" 2 . 

- ionized reflection . We replaced the Compton-reflection 
component in M pr with a model explicitly including the 
ionization state of the refiectoi0, producing also the ion- 
ized iron lines which in other models are fit empirically 
with Gaussian profiles. This model is based on the calcula- 
tion by Ross & Fabian (2005). If we assume rj„t=1.8, the 
model does not yield a good fit to the best-quality spectra 
(see, e.g., the results on the deepest XMM-Newton observa- 
tion in Fig. m . If left free to vary, the spectral index pegs to 
the minimum value allowed by the model (1.6), still leaving 
significant residuals. Even when applied on spectra with a 
lower statistical quality, where the formal quality of the fit 
is comparable to that yielded by M pr , the model exhibits: 
a) too broad an iron line profile, and; b) too shall a photo- 
electric absorption edge with respect to the data. We will 
not consider therefore this model any longer in this paper. 

In Fig. [5] we compare the light curves of the 4-5 keV flux 
due solely to the Compton-reflection component as mea- 
sured in different spectral models, and corrected for the 
intervening absorption in the M or model. Differences are 
apparent. The most evident is that the error bars in M or 
and M urot are larger, due to the larger number of degrees of 
freedom. However, the overall Compton-reflection variabil- 
ity pattern is not qualitatively modified. The main change 
is the measured flux of the March 2012 EPIC-pn observa- 
tion. It is larger by about 0.1 dex with M or when compared 
to M pr , and is fully unconstrained with M urot . 



5 through the ionization parameter £ = L/(nR 2 ), where L is 
the ionizing luminosity, n is the electron density, and R the dis- 
tance between the reflector and the source of ionizing continuum 
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Fig. 5. Light-curve of the absorption-corrected 4-5 keV 
flux of the Compton-reflection component when models 
M pr (upper panel; the same quantity as in the upper panel of 
Fig- [Hi Mor {medium panel), and M urot (lower panel) are 
used. The units on the y-axis are the difference between 
the logarithm of the flux and its median. The dotted line 
indicates the median level. 



In summary: the fact that the flux associated to the 
Compton-reflection component varies does not depend on 
the spectral model used to fit the data analyzed in this 
paper (within the simple spectral models used). On the 
other hands, quantitative estimates of the variability pat- 
tern timescales do depend on the spectral modeling. Reader 
should bear in mind this important caveat over the rest of 
the paper. Long-term deeper monitoring of Markarian 3 
would be required to break this degeneracy, where each in- 
dividual spectrum has got a statistical quality comparable 
or better than the deepest EPIC-pn observation discussed 
in this paper. 

Finally, in the lower panel of Fig. [1] we show the light 
curve in the 0.5-1.5 keV energy band. It was measured by 
fitting the spectra in the 0.5-2 keV with a phenomeno- 
logical model constituted by a power-law (of photon index 
r s ) and 7 unresolved Gaussian profiles, with energy and 
normalization fixed to the best-fit values obtained from 
the deepest XMM-Newton observations. We did not in- 
clude the Suzaku measurement in this panel, because it is 
likely to be contaminated by the ROSAT source IXO30 
(jBianchi et al. 20051 lAwaki et al. 2008[) . Variability with 
a dynamical range of about ±10% is seen also in this 
band (Xg/V=249.0/12). Such a dynamical range is of the 
same order of the systematic uncertainties in the cross- 
calibration among flying X-ray detectors in this energy 
range ( |Tsujimoto et al. 201 ip . This casts doubts on an as- 
trophysical origin of these differences. We discuss this find- 
ing in more details in Sect. 15.41 of this paper. 



4. Variability timescales in Markarian 3 

In Fig. [5] we compare the 4-5 keV observed flux against the 
flux measured by the Swift/BAT in the 14-195 keV energy 
band. Both curves are normalized to their respective mean. 
The 4-5 keV points are calculated with the M pr model. 

The variability observed in the BAT light curve could 
be due either to true changes of the intrinsic nuclear flux 
or to changes of the column densities covering the source of 
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Fig. 7. Swift/BAT Markarian 3 light curves in three bands, 
and their hardness ratio. 



the primary emission. In the latter case we should see differ- 
ent variability patterns in light curves extracted in different 
BAT energy ranges, with the light curve extracted from the 
hardest band tracing most closely the primary continuum. 
In Fig. [7]we show the Markarian 3 BAT light curves in three 
bands: 15-30 keV; 30-70 keV, and 70-150 keV, together 
with their hardness ratios. They exhibit a qualitatively sim- 
ilar pattern. The hardness ratios are consistent with being 
constant, with x^/v values of 19.7/43 and 22.3/41 for the 
30-70/15-30 and the 70-150/30-70 keV hardness ratio, re- 
spectively. 

With the light curves shown in Fig. [6] and [7] we can 

estimate the following timescales: 

- The minimum variability timescale of the primary emis- 
sion, t p . We define it from the shortest interval be- 
tween two bins of the BAT light curve which are in- 
consistent at the 95% confidence level. This quantity 
is basically constrained by the statistics of the BAT 
light curve. There are several consecutive bins, both in 
the full as well as in the energy-resolved BAT light 
curves, which are inconsistent at this statistical level. 
We conclude that r p <515 days. Significantly shorter 
(hours to days) timescales are common in Seyfert galaxy 
dBarr k Mushotzky 19861 |McHardy et al. 2006D , but they 
are basically inaccessible to us in such heavily obscured 
source as Markarian 3 

- The maximum variability timescale of the Compton flux, 
t. We define it as the minimum time interval between two 
bins of the Compton flux light curve differing at more than 
the 95% confidence level. It depends on the spectral model, 
due to the larger statistical errors affecting the measure- 
ment of the Compton flux where more complex spectral 
models are used (Tab. [3]) 

- The delay in the response of the Compton flux to variation 
of the primary continuum At s We calculate the quantity: 



X 2 (At s ) = E[(F 4 "_ 5 (Ai s ) - F^_ 195 )/(*l_ 5 + o? 4 _ iaB )] 

where: F£ E are the normalized fluxes in the AE energy 
range; ctab are the corresponding statistical uncertainties, 
augmented by 5% to take into account possible systematic 
cross-calibration uncertainties; and At s is a temporal shift 
applied to the 4-5 keV flux data points only. For this cross- 
correlation to be astrophysically meaningful (cf. Sect. 15.1)) . 
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Fig. 6. Light curve in the 4-5 keV energy band (filled circles) against the Swift/BAT light curve in the 15-150 keV 
energy band (empty squares). Both curves are normalized to their respective mean for direct comparison. 
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Fig. 8. Probability of x 2 (At s ) as a function of At s for M pr 
when AtBAT = 75 days. The dotted line indicates the 95% 
confidence level. 



Table 3. Minimum variability timescale of the Compton- 
reflection component, r, and delay in the response of the 
Compton-reflection to changes of the primary continuum, 
At s . All quantities are in days. Empty cells indicate uncon- 
strained measurements. 



Model 


r 


At s 


M pr 


64.0 ± 0.3 a 


> 1162.5 


M or 


152.2 ±0.2 b 






115.64 ±0.13 c 





a time difference between the 2012-01-11 and the 2012-03-15 ob- 
servations. 

fc time difference between the 2000-10-19 and the 2001-03-20 ob- 
servations. 

c time difference between the 2012-01-11 and the 2012-06-05 ob- 
servations. 



the BAT light curve must be binned with At bat > t. We 
used the smaller integer multiple of 15 days fulfilling this 
relation. 

We consider this cross-correlation statistically signifi- 
cant when the one-sided probability of x 2 (Ai s ) is larger 
then 95%. For M pr , a family of solutions exists, correspond- 
ing to At s £ 1200 days (see Fig.©. If the 4-5 keV Compton 
flux is calculated using the other spectral models, no At s 
satisfies our statistical criterion. 

In principle, the variable Compton-reflection compo- 
nent could contribute to the variability of the BAT light 
curves as well. While spectral decompositions in contin- 
uum components must been taken with a grain of salt 
when performed in spectral regions where the instru- 
mental resolution is poor, the results of the Suzaku ob- 



servation formally indicate that the Compton-reflection 
could contribute about one third of the total flux at 
15 keV, and about 5 0% thereof above 30 keV (cf. Fig.2 
in lAwaki et al. 2008)) . However, it is unlikely that the 
Compton-reflection contribute significantly to the variabil- 
ity above 15 keV on the following grounds: a) the short 
variability timescale of the BAT light curve; b) its energy 
invariance: Compton-reflection and primary continuum ex- 
hibit very different spectral shape around the Compton 
peak at 30 keV; c) the fact that there is very little correla- 
tion between the 4-5 keV and the 15-150 keV light curve at 
zero time shift. The probability of x 2 (At s = 0) is ^6xl0~ 4 
for M„ r . 
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5. Discussion 

5.1. Constraining the torus geometry through variability 

In principle the variability timescales determined in Sect.0] 
can be used to constrain the location and geometry of the 
optically-thick reprocessor. We discuss in this Section the 
method, its assumptions, and its limitations, as well as the 
results of its application to the monitoring campaign of 
Markarian 3 discussed in this paper. 

If the primary continuum illuminating a slab of optically 
thick gas varies on a timescale t p , the variability of the 
primary emission can be fully diluted if the light crossing 
time across the reflecting part of the slab is larger than t p . 
For all the spectral model discussed in this paper r 3> r p 
(Tab. [3]). Due to its finite size, the reprocessor will not 
respond to the instantaneous primary flux, rather on an 
average over its own light-crossing time. This is the reason 
why we impose in our cross-correlation analysis (Sect. [4]) 
the condition AtsAT > T. 

Under this condition, the delayed response of the re- 
flected component to changes of the primary continuum 
can be used to constrain the distance between the AGN 
and the reflecting cloud, through simple light crossing-time 
arguments: Al ~ cAt Sl where Al is the difference in opti- 
cal path between the primary and the reflected radiation. 
If reflection occurs from the inner far side of the torus, 
d ~ AZ/2, where d is the distance between the reflecting 
cloud and the source of illuminating continuum. 

If reflection occurs in an individual cloud the difference 
in optical path length between photons coming from ex- 
treme regions of the cloud has to be smaller than its light- 
crossing cloud linear size. This puts a lower limit on its 
opening angle: 



> arcsin( 



1 



1 + S/d 



where S is the size of the reflecting cloud and d its distance 
from the primary source. For M pr : S <0.05 pc, d <;0.5 pc, 
0^71°. 

An independent estimate of the reflector opening an- 
gle can be derived from the opening angle of the ion- 
ization cone, 9 . Capetti et al. (1995) derive 8 >55° 
from HST imaging of the NLR. This estimate is larger- 
then previous estimates based on ground-based observa- 
tions (| Wilson k Tsvetanov 1994|) . 

Ikeda et al. (2009) tried to constrain the geometry of 
the reflector by fitting the Suzaku time-average spectrum of 
Markarian 3 with Monte-Carlo simulations of the primary 
continuum emission by a Compton-thick accreting torus. 
The simultaneous measure of the reflected and of the trans- 
mitted (above 10 keV) components may in principle allow 
to constrain the geometry as well as the column density 
of the reflector ( |Murphy k Yaqoob 2009) |Yaqoob 2012 1. 
Unfortunately, there is an intrinsic degeneracy between the 
opening angle and the inclination angle against the line-of- 
sight, which cannot be resolved even with the best-quality 
broad band spectral data available. Awaki et al. (2008) con- 
strained the distance of the reflector within ~1 pc due to the 
low- ionization and unresolved profile of the iron line. A sim- 
ilar analysis of the deepest EPIC-pn spectrum by Bianchi 
et al. (2005) yielded an estimate of the inner size of the 
torus comprised between 0.3 and 0.5 pc. 

While the formal geometrical constrains on the torus 
derived by our variability study are not inconsistent with 



those obtained by other studies using different techniques, 
one might ask whether the implicit underlying picture of 
a single compact reflecting cloud with a size smaller them 
0.05 pc is tenable. We present in the following further ob- 
servational evidence which questions it. 



5.2. The K a iron line emission region 

In the current picture of Compton-thick reprocessing in 
the AGN environment, the Compton- reflection continuum 
and the "narrow"0 iron K Q fluorescent line should be pro- 
duced by the same reflecting gas. Indeed, these features 
almost invariably come together in the X-ray spectra of 
highly-obscured AGN. Existing exceptions (e.g., NGC7213, 
Bianchi et al. 2008) are explained by reprocessing from 
Compton-thin gas in BLR clouds. If both components are 
produced by reprocessing of the same emission, one may 
expect them to respond simultaneously to changes of the 
primary. 

This does not seem to be the outcome of the X-ray 
spectroscopic campaign in Markarian 3. The intensity of 
the strong and ubiquitous iron K a fluorescent line does 
not follow the variability of the Compton-reflected contin- 
uum (Fig. [5]). This finding is puzzling. In order to explain 
it, one has to assume that the production of the line is 
decoupled from the underlying continuum. The line could 
be produced in Compton-thin clouds, while the Compton- 
reflection continuum could track Compton-thick clouds. 
Several mechanisms can be invoked to explain the lack of 
a strong iron emission line from Compton-thick material. 
Ionization can suppress the emission line via resonant trap- 
ping (IMatt et al. 1 993) or full electron stripping. However, 
the data with the best statistical quality rule out this expla- 
nation (cf. Fig.@|. Alternatively, the Compton-thick clouds 
could be iron deficient. It is, however, unclear why abun- 
dances of different gaseous systems in the environment of 
the same AGN should be largely different. Furthermore, 
this hypothesis is inconsistent with the deep photoelec- 
tric absorption edge, clearly visible in all spectra (see, e.g., 
Fig. ITU]). More importantly, there is no other evidence of the 
presence of Compton-thin clouds from X-ray spectroscopic 
observations of Markarian 3, although they could not be 
detected along the same line-of-sight of the Compton-thick 
absorber. 

To further investigate this point, we have reanalyzed 
the 0*^ order image of the Chandra /HETG observation of 
Markarian 3 originally presented in Sako et al. (2000). As 
they showed (cf. the Fig.l in their paper), the soft X-ray 
emission is extended along an approximate E-W direction 
(the elongation direction corresponds to an astronomical 
position angle of 356°. 6). Intensity profiles in three differ- 
ent bands: 0.5-3 keV, 3-6 keV, and 6-7 keV are shown in 
Fig. |51 We confirm that emission is primarily extended in 
the softest energy band along the E-W direction up to ~4" 
on the East side and ~2" on the West side. In the hard band 
(E > 3 keV) an excess by ±1" (~270 pc) is visible in the 
hardest bands along the North direction only (i.e., perpen- 
dicular to the extension of the X-ray NLR). However, there 



6 We use the term "narrow" here to identify emission line fea- 
tures, whose profile is not broadened and skewed by relativistic 
effects in X-ray illuminated accretion disks a few gravitational 
radii from the black hole event horizon (Fabian et al. 1989|) 
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Fig. 9. Intensity profiles measured from the t/l order 
image taken during the Chandra/HETG observation of 
Markarian 3 (jSako et al. 2000|) . From top: E-W; bottom: N- 
S, Symbols indicate energy bands: empty circles: 0.5-3 keV; 
filled squares: 3-6 keV; filled triangles: 6-7 keV. The dotted 
lines indicate the Point Spread Function. From the inner- 
most to the outermost profiles in each panel: 0.5-3 keV, 
3-6 keV, and 6-7 keV, respectively. 

is no statistical difference between the Compton continuum 
and the iron fluorescence line profiles. 

Extended soft X-ray (E < 2 keV) emission on scales 
ranging from a few hundreds parsecs to a few thousands 
parsec has been commonly discovered in highly obscured 
nearby AGN, almost invariably in morphological agreement 
with the optical NLR ( Bianch Tet al. 20061 and references 
therein). Extension in the hard energy band is much less 
common. Young et al. (2001) reported hard X-ray extended 
emission, together with a tentatively detected iron line, in 
NGC 1068 on scales ~20" (2.2 kpc), alongside a low sur- 
face brightness component up to ~50" (5.5 kpc). More re- 
cently, Marinucci et al. (2012b) reported a clumpy and flat- 
tened ^150 pc structure in the hard band Chandra/ ACIS 
image of NGC 4945. Compact hard X-ray emission in 
the Circinus Galaxy (jSambruna et al. 200ip . constrains the 
size of the Compton reflector to within 8 pc. In bright 
nearby obscured AGN, mir-IR intcrferometry allows to con- 
straint the location and size of the dusty torus within 
a few parsecs (jJaffe et al. 2004( [Meisenhoim er et al. 2007, 
|Tristram et al., 2007} . While the X-ray and IR studies in- 
dicate that dusty structures in the nuclear environment of 
AGN are extended on scales from a fraction to hundred par- 
secs, we stress that the two techniques probe two different 
phases of the obscuring matter in the nuclear environment. 
X-ray direct imaging and variability are sensitive to the 
gaseous phase, whereas IR measurements probe the dust 
component. 

5.3. Implications for the AGN structure model 

The results presented in this paper favour a com- 
plex structure of the optically-thick absorber/reflector in 
Markarian 3. X-ray energy-dependent cross-correlation re- 
sults favor a reflecting cloud at a distance <;0.5 pc. On the 
other hand, the variability of the Compton-reflection com- 
ponent on scales ^64-150 days (depending on the spectral 
model) suggests that the variability episodes presented in 



this paper might be dominated by small clouds, implying 
a large opening angle, i.e. a small covering factor. We dis- 
covered a puzzling lack correlation between the Compton- 
reflection flux and the Fe K Q line intensity from our X- 
ray monitoring (Fig. [3J correlation probability <510 -4 ). 
Furthermore, at least part of the reflection occurs on scales 
as large as 300 pc (Fig. . Indirect evidence for optically 
thick reprocessing of the primary radiation with a large 
covering fraction comes from the fact that Markarian 3 
is a strong mid-IR emitter. It is one of the prototypi- 
cal AGN-dominated Type 2 objects in mid-IR continuum 
(|Deo et al. 2009|) . This also represent an indirect evidence 
in favour of an extended structure reprocessing the AGN 
emission into the IR band. 

In summary, we report in this paper various (and partly 
conflicting) evidence for optically-thick reprocessing of the 
AGN primary emission in Markarian 3 on a variety of dif- 
ferent spatial scales, from fraction to hundreds of parsecs. 

Comparing these findings with other Compton-thick 
AGN in the local Universe unveils a similar variety, 
although detecting variability of the Compton-reflection 
component is difficult and therefore rare. In NGC 4945, 
the historical (Chandra, Suzaku and XMM-Newton) 
light curve of the reflected component remains constant 
within <10% despite strong intrinsic variability of the 
Swift/BAT light curve on time scales longer than 1 year 
(jMarinucci et al. 2012b[) . Constant flux of the reflected 
components over a monitoring campaign covering almost a 
decade was reported in NGC 5506 (Guaina zzi et al. 2010. 
On the other hand, a delayed response of the K a iron 
line and Compton-reflection to a decrease of the pri- 
mary flux by a factor ~20 in 16 years constrained the 
torus to be located at about 3.2 pc from the AGN in 
NGC 2992 (jWeaver et al. 1996|) . Variability by about 25% 
of the iron K a line on timescales of about 1 years were 
reported by different authors in the unobscured Seyfert 
NGC 4151 UPerola et al. 19861 ITakahashi et al. 20021 
Zdziarski et al. 2002, Schu rch et al. 2003)) . suggesting that 
the emitting region is at about 0.03 pc from the nucleus. 

One possible explanation of these apparently contra- 
dicting results is that the absorber is clumpy. Qualitatively, 
reflection on different spatial scales require that absorption- 
free line-of-sights to the active nucleus must exist for clouds 
located at various distances. The case for a clumpy torus 
was originally advocated by Krolik & Begelman (1998), 
and more recently fully developed in the calculations by 
Nenkova et al. (2002, 2008). While strong direct observa- 
tional evidence exists for a clumpy absorber on scales of the 
BLR in several nearby AGN from X-ray variability studies 
([Elvis et al. 20041 IRisaliti et al. 20051 IPuccetti et al. 20071 
Bianchi et al. 2009b, IMaiolino et al. 20101 

Risaliti et al. 2011), evidence for the torus dumpiness 
are mainly yielded by the comparison between torus mod- 
els and multiwavelength Spectral Energy Distributions. 
A compact, pc-scale torus can h ardly reproduce th e 
multi-temperature IR emission (|Pier &: Kro lik 1992). 
Large-scale dusty torii can reproduce the observations 
(|Granato &: Danese 1994p . but are nowadays disfavoured 
given the IR interferometric evidence of very compact 
warm dust structures sur rounding nearby obscur ed AGN 
(see, however, NGC 4945: IMarinucci et al. 2012bp . 

Recently, Eliztur (2012) proposed that the complex 
structure of the absorber/reflecting matter in the AGN 
environment should lead us to a revision of our interpre- 
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tation of the AGN Unification Scenario. Its strictest for- 
mulation postulates that the orientation with respect to an 
azimuthally-symmetric obscuring structure is the only driv- 
ing factor in the classification of an AGN. Elitzur suggests 
that this view should be replaced by a probabilistic inter- 
pretation, whereby the classification of a radio-quiet AGN 
would depend on the covering factor distribution of the ob- 
scuring matter. "Type 2" would be more likely identified in 
AGN surrounded by matter with a higher covering factor. 
A least extreme view, where the Poissonian fluctuations in 
the number of optically-thick clouds on scales comparable 
or lower then the BLR is superposed to the standard formu- 
lation of the Unified Scenario is discussed by Matt (2000) 
and Bianchi et al. (2009). 

Markarian 3 is one of the few obscured AGN in the lo- 
cal universe, where evidence of AGN primary emission re- 
processing by optically thick gas on scales from a fraction 
of to several hundreds parsecs have been reported. It may 
represent one of the best laboratories to test these scenar- 
ios with future observations in X-rays as well as in other 
wavelengths. ALMA is expected to give a paramount con- 
tribution to our understanding of the circumnuclear dusty 
torus. 



5.4. On the variability in the soft X-ray band 

The soft X-ray [E < 3 keV) flux history in Markarian 3 
also exhibit significant variations around a mean of 
~7xl0 -13 erg cm" 2 cm" 1 . The dynamical range is 
around ±10%. This is of the same order of magnitude 
as the typical systematic uncertainties in the flux deter- 
mination between the EPIC-pn and the XRT detector 
flTsujimoto et al. 2011[ ). Indeed, the variability pattern is 
dominated by the comparison between the EPIC-pn and 
the XRT observations, although the difference (XRT yield- 
ing systematically lower fluxes) is opposite to what mea- 
sured by |Tsujimoto et al. 201 1| A 5% difference in the X- 
ray flux measured by XMM-Newton on 2000-2002 and 2012 
could be affected by uncertainties in the background sub- 
traction in the most recent observation (see Sect. 2). While 
we stress that further higher-quality measurements are re- 
quired to confirm this finding, we discuss briefly in this 
Section some of its possible implications. 

Bianchi et al. (2006) estimate that about 10% of the 
emission in this band is unresolved in the Chandra/ ACIS 
image. A variability by a factor ~100% of this com- 
pact component (associated to the innermost region of 
the NLR) could account for the observed variability pat- 
tern. Alternatively, the observed variability could be due 
to the spilling of the primary emission through a par- 
tial covering absorber. In the Palermo BAT Catalogue 
(jCusumano et al. 2010|) Markarian 3 is detected with a 15- 
150 keV flux of (1.18 ± 0.03) x 10" 11 erg cm" 2 cm" 1 and a 
spectral index T bat = 1-74. The extrapolation of this flux 
into the 0.5-1.5 keV energy band, assuming a simple power- 
law of photon index Vbat is ~4.4xl0~ n erg cm" 2 cm" 1 . 
The observed flux is about 1.5% of this extrapolation. 
Such a fraction is not unusual in Compton-thin AGN 
(jRisaliti 2002|) . See also Paggi et al. (2012) for the detection 
of variability in the soft X-ray spectrum of Mkn 573. 



6. Conclusions 

The main conclusions of this paper can be summarised as 
follows: 

— we report on a study of the decade-long X-ray flux vari- 
ability in Markarian 3 in the 4-10 keV energy band, 
dominated by Compton-reflection of the primary AGN 
emission. The shortest variability timescale is ^64 days, 
as measured by the shortest interval between two sta- 
tistically inconsistent measurements 

— assuming that the light curve of the Compton-reflection 
dominated energy band and that of the 15-150 keV en- 
ergy band are correlated, with the former lagging the 
latter, the minimum measured delay is *J1200 days 

— the intensity of the iron K a emission line is uncor- 
rected to the flux in the Compton-reflection domi- 
nated band. This suggest that the bulk of the iron 
line is emitted in a different, more extended structure. 
Energy-resolved analysis of existing Chandra data indi- 
cates a region of extended emission of both the line-free 
Compton-reflection continuum and of the Fe K Q line up 
to ^300 pc to the North of the active nucleus 

Markarian 3 exhibit a rich phenomenology associated 
to the optically-thick reprocessing of the AGN primary ra- 
diation. Evidence of reprocessing on very different spatial 
scales have been presented and discussed in this paper. 
They flow with the stream of mounting experimental ev- 
idence in favor of an extended and clumpy structure of the 
absorbing/reprocessing gas and dust in the innermost par- 
secs around the super-massive black hole. Revisions of the 
standard interpretation of the Unified Scenario encompass- 
ing this com plexity have been suggest ed by several authors 
(jMatt 2000UElvis 20001 [Elitzur 2012|) . 

Still, the constraints on the geometry of the reprocess- 
ing matter in Markarian 3 are dependent on the details of 
the spectral modeling, due to the coarse X-ray monitoring 
in the Compton-reflection dominated energy band, as well 
as to the statistical quality of the available data. We intend 
to tighten them in the nearby future, by pursuing a specific 
observational monitoring program. The perspectives offered 
by the successful launch of NuSTAR (jHarrison et al. 2010]> 
are particularly promising. It is reasonable to assume that 
at least part of the material reflecting the intrinsic contin- 
uum is also responsible for the absorption along the line-of- 
sight. In this case variability of the measured line-of-sight 
may provide additional constraints. These measurement are 
virtually impossible now, but may be possible thanks to the 
increased sensitivity of NuSTAR at energies larger then the 
photoelectric cut-off. Similar measurements may be possi- 
ble also with Astro-H in a nearby future. 

7. Appendix A 

In Fig. [TO] we show the spectra {upper panels) and residu- 
als in units of data/model ratio (lower panels) for all the 
observations analyzed in this paper. 
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Fig. 10. Spectra {upper panels) and residuals in units of data/model ratio {lower panels) when the baseline model is 
applied to the spectra of the observations discussed in this paper. The data were rebinned such that each data point 
corresponds to a sigma-to-noise ratio >3 for visualization purposes only (the analysis was performed on the unbinned 
spectra). The spectrum of the deepest observation (October 2000; upper left) is shown with a coarser binning (minimum 
signal-to- noise ratio > 10) in Fig. [4j 
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